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Abstract 

Background: Transcription factors fronn tlie MADS-box family play a relevant role in cell differentiation and 
developnnent and include the aninnal SRF (serunn response factor) and 1\/1EF2 (myocyte enhancer factor 2) proteins. 
The social amoeba Dictyostelium discoideum contains four genes coding for MADS-box transcription factors, two of 
these genes code for proteins that are more similar to SRF, and the other two code for proteins that are more 
similar to MEF2 animal factors. 

Results: The biological function of one of the two genes that codes for MEF2-related proteins, a gene known 
as mef2A, is described in this article. This gene is expressed under the transcriptional control of two alternative 
promoters in growing cells, and its expression is induced during development in prespore cells. Mutant strains 
where the mef2A gene has been partially deleted were generated to study its biological function. The mutant 
strains showed reduced growth when feeding on bacteria and were able to develop and form fruiting bodies, 
but spore production was significantly reduced. A study of developmental markers showed that prespore cells 
differentiation was impaired in the mutant strains. When mutant and wild-type cells were set to develop in 
chimeras, mutant spores were underrepresented in the fruiting bodies. The mutant cells were also unable to form 
spores in vitro. In addition, mutant cells also showed a poor contribution to the formation of the tip-organizer and 
the upper region of slugs and culminant structures. In agreement with these observations, a comparison of the 
genes transcribed by mutant and wild-type strains during development indicated that prestalk gene expression was 
enhanced, while prespore gene expression decreased in the mefIA strain. 

Conclusions: Our data shows that mefIA plays a role in cell differentiation in D. discoideum and modulates the 
expression of prespore and prestalk genes. 
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Background 

Mef-2-related transcription factors belong to a family of 
proteins that are present in all eukaryotic organisms [1,2]. 
These proteins share a very conserved DNA-binding and 
protein-dimerization domain, the MADS-box, named 
after the transcription factors MCMl (from yeast), 
Agamous, Deficiens (from plants) and SRF (from animals) 
[3,4]. Two subfamilies of MADS-box transcription factors 
have been defined according to the MADS-box sequence: 
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type I and type II [5]. Plants have a large number of types 

I and II MADS-box proteins while other organisms, such 
as fungi and animals, usually have one or more proteins of 
each subfamily [6]. Animals, for example, have only one 
type I protein (serum response factor [SRF]) and four type 

II proteins (myocyte enhancer factor 2 A-D [Mef2 A-D]). 
The two types of factors recognize different A/T rich bind- 
ing sites. SRF and related factors recognize the consensus 
sequence CC{KIT)^GG [7,8], while Mef-2-related factors 
recognize the CTA(A/T)4TAG consensus sequence [9]. 

MADS-box transcription factors accomplish various 
biological functions. In plants, they are critically involved 
in floral formation and development [5]. In yeasts. 
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MCMl participates in the regulation of pheromone 
expression, metabolism [10] and DNA replication [11]. 
In animals, these transcription factors are mainly 
involved in the regulation of cell- differentiation processes. 
SRF deletion is lethal in mice because cells are impaired in 
cell adhesion and migration [12], and the embryo cannot 
complete gastrulation [13]. Tissue-specific deletion 
has shown that SRF is also required for terminal dif- 
ferentiation of skeletal, cardiac and smooth muscle 
cells and for neural cell migration [14]. Additional studies 
have demonstrated that SRF regulates the expression of a 
large number of genes coding for actin-cytoskeleton-related 
proteins [15,16]. 

Mef-2 proteins are involved in regulating the expression 
of muscle-specific genes, both in Drosophila [17] and in 
mammals [18], in collaboration with MyoD-related tran- 
scription factors [19]. In mammals, there are four genes 
coding for very similar Mef-2 factors that appear to 
functionally complement each other, at least partially 
(Mef2A-D) [20]. However, Mef2C-null mice die early in 
their development due to cardiovascular abnormalities [21], 
and Mef2A-null mice die perinatally from heart defects 
[22]. In addition, numerous studies have shown that Mef-2 
factors are also involved in the differentiation of several 
other cell types, such as neural crest cells, endothelial cells, 
chondrocytes, neurons and lymphocytes [23,24]. 

Our group has approached the functional study of 
MADS-box transcription factors in the social amoeba 
Dictyostelium discoideum. These unicellular organisms 
live in forest soils, feeding on bacteria and other microor- 
ganisms and are able to develop as multicellular organ- 
isms under starvation conditions. In such conditions, up 
to 10^ individual amoebas aggregate to form a fruiting 
body composed of a basal disk, stalk and sorus where up 
to 80% of the original amoeba differentiate into resistant 
forms called spores [25-27] . The initial step is the aggrega- 
tion of the cells towards cAMP-secreting centers to form 
a mound. Cells within the aggregates initiate a differenti- 
ation process to form two main cell types: prestalk and 
prespore. Prestalk cells migrate to the top of the mound, 
emerging as a tip. A culmination process is later initiated 
by the migration of prestalk cells from the tip towards the 
substrate through the mass of prespore cells, piling up 
and terminally differentiating to form the stalk. The mass 
of prespore cells remains attached to the top of the 
forming stalk, rising from the substrate until culmination 
is completed. Migratory structures, called slugs, can be 
formed before culmination under adverse environmental 
conditions. In this case, the slugs migrate towards warmer 
and lighter places for culmination to facilitate the dissemi- 
nation of the spores. By the end of culmination, prespore 
cells differentiate inside the sorus to form mature spores. 

Analysis of the D, discoideum genome has shown that it 
contains four genes coding for MADS-box transcription 



factors, namely srfA, B, C and D. Previous studies in our 
laboratory have shown that srfA is required for the proper 
development of the fruiting body, including the slug 
migration and culmination steps, and is essential for spore 
terminal differentiation [28,29]. The srjB gene is expressed 
earlier than srfA during development, and the encoded 
protein is involved in the initiation of the developmental 
process, cell migration and the initiation of culmination 
[30]. The functional study of the srfC gene is described in 
this article. We present evidence demonstrating that srfC 
is more similar to animal Mef-2 genes and propose 
naming it mefIA, By analyzing the phenotype of mutant 
strains and gene expression levels during development, we 
clearly show that this protein is involved in D, discoideum 
development and, in particular, in the differentiation of 
prespore cells and one group of prestalk cells. 

Results 

Characterization of the meflJK (srfQ gene 

The analysis of the D, discoideum genome identified four 
genes coding for proteins with regions similar to the 
MADS-box domain. These genes were named srfA, B, C 
and D. The putative MADS-box region of the proteins 
encoded by these genes was analyzed in more detail and 
compared to that of vertebrates (//. sapiens^ G. gallus, X. 
laevis), invertebrates (D. melanogaster, A. franciscana), 
amoeba (£. histolytica) and fungi {S. cerevisiae) type I (SRF, 
MCMl, ARGl) and type II (Mef2) proteins. Amino acid 
sequences were aligned using the ClustalW and ClustalX 
programs (Additional file 1: Figure SI). The multiple align- 
ments were used to calculate the phylogenetic tree shown 
in Figure 1. The results indicate that the D, discoideum 
genes srfA and srfB code for proteins that are more similar 
to type I genes, such as the animal SRF gene, than to type 
II genes, such as Mef2. D, discoideum SrfA and SrfB appear 
to form a monophyletic group that more closely resembles 
to animal proteins than to fungi proteins. In contrast, the 
proteins encoded by D, discoideum srfC and srfD are more 
similar to animal MEF2 than to SRF proteins. In this case, 
the protein encoded by srfC is more closely related to that 
of the amoeba E, histolytica and is also more similar to 
animal Mef2 proteins than to SrfD. In fact, SrfD appears to 
have diverged significantly from the other type II proteins 
analyzed. The results of this analysis centered our attention 
on srfCy which we propose to rename as mefZA, a name 
that we will use throughout the rest of the article. 

The expression of mef2A during vegetative growth 
and development was analyzed by RT-PCR, and the 
results are shown in Figure 2A. Expression was 
detected in growing cells (time 0 in Figure 2A), but a 
large induction was observed at 4 hours of development, 
an induction that was maintained at later developmental 
stages. The promoter region of the gene was characterized 
to further determine the temporal and spatial patterns of 
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Figure 1 Phylogenetic tree of the D. discoideum MADS-box-containing proteins. The amino acid sequences of tine MADS-box region from 
MEF2 and SRF proteins from vertebrates {H. sapiens, X. loevis, G. gollus), invertebrate (D. melonogoster, A. fronciscono) animals, fungi (S. cerevisioe) 
and amoeba {E. histolytica) were compared to tliose of tine four D. discoideum MADS-box containing proteins (SrfA, B, C and D) using tine ClustalW 
program at tine online Biology Work Bench facility from the San Diego Supercomputer Center (http://workbench.sdsc.edu). Phylogenetic trees were 
determined using the neighbor-joining method and the ClustalX program. A random generator seed of 1 1 1 was used, and 1000 bootstrap trials were 
calculated. The number of times that each branch was obtained is indicated at the base of each branch. The tree was drawn using the NJplot 
program. The evolutionary distance scale, calculated as a fraction of amino acid changes, is indicated at the left lower corner of the figure. 



expression. Primer-extension experiments were performed 
using the RACE (rapid amplification of cDNA ends) tech- 
nique. The data indicated the existence of one intron in 
the 5' UTR and two regions of transcription initiation 
(Additional file 1: Figure S2) [GenBank:KC852901]. The 
gene structure diagram, which includes the 5' untranslated 
region (UTR) and the transcription initiation sites, is shown 
in Figure 2B. The more upstream transcription initiation 
sites of each region are located at nucleotides -144 
and -814, respectively. The first intron is located in 
the 5' UTR of one of the two transcribed mRNAs 
(nucleotides -536 to -371, in relation to the A of the 
initiation codon). The data obtained from the experiments 
also enabled us to define the limits of the second intron, 
located between the first (A) and second nucleotide (T) of 
the translation initiation codon. 

The transcriptional activity of both promoters was ana- 
lyzed by the use of reporter vectors where Promoter 1 
(Prl), Promoter 2 (Pr2) or the complete promoter region 
(cPr) were cloned, thereby driving lacZ expression. Pools 
of transformed cells obtained for each promoter were ana- 
lyzed for p-galactosidase activity. Promoter 1 drove lacZ 
expression in scattered cells at the mound and finger 



stages of development, but its activity markedly increased 
in the prespore region of slug, Mexican-hat and culminant 
structures (Figure 2C). Promoter 2 was active in scattered 
cells of aggregates and fingers, but the activity decreased 
almost completely at later developmental stages, except 
for a few cells in the basal disk of culminant structures. 
The activity of the complete promoter showed the sum of 
Prl and Pr2 and was maximal in the prespore region of 
developing structures. 

Generation of mef2>A-deficient strains 

The study of the biological function of mefZA was 
approached through the generation of mutant strains 
where the gene was partially deleted by homologous re- 
combination. The deleted region included the first two 
exons, coding for the S' untranslated region of the gene, 
and the third exon, coding for the N-terminal region of 
the protein, including the MADS-box domain (Figure 2B). 
Several mutant clones were isolated on two different back- 
grounds, the AX2 and AX4 D, discoideum axenic strains. 
Figure 3A shows the RT-PCR analyses demonstrating 
mep,A gene deletion in one AX4-derived clone (37) and 
two AX2-derived clones (2, 3). Later results show no 
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Figure 2 Structure and expression of the mef2A gene. Panel A. RNAs were isolated from AX4 cells during growth (0) or after development 
(2-24 hours). Expression of mef2A or the large mitochondrial rRNA (internal control), were determined by PCR. Panel B. RNA purified from AX4 
cells developed for 8 hours was converted into cDNA, and the 5' end extended using an oligonucleotide complementary to nucleotides 145 to 
164 of the gene. Two different amplification products were obtained and their nucleotide sequences aligned to the genomic DNA. A diagram of 
the deduced structure of the gene is shown. Exon sequences are indicated as boxes, blue boxes for untranslated regions and deep red boxes for 
translated regions. Arrows indicate transcription initiation sites. The sequence has been numbered from the A of the translation initiation codon. 
The location of intron/exon borders is indicated on the lower part of the diagram and that of transcription initiation sites in the upper part. The 
DNA fragment deleted in the mutant strains is indicated in the lower part of the diagram. Panel C. The activity of the mef2A promoter regions 
was studied using locZ reporter vectors driving the expression of a short-lived form of (3-galactosidase. The region from the 3' end of the closest 
upstream gene (-2201) to the end of exon 1 (-489) (Prl) and from the end of exon 1 (-51 1) to exon 3 (164) (Pr2), or the complete promoter 
region (-2201 to 164) (cPr), were cloned. /acZ expression was analyzed by histochemistry, using the Xgal substrate (2 hours of incubation) in 
aggregates (10 hours of development), finger (16 hours), slug (24 hours of development under migration conditions), Mexican hat (18 hours) and 
culminant (22 hours) structures. Pictures were taken using a Leica stereomicroscope, after counter-staining with eosin. Scale bar: 0.2 mm. 
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Figure 3 Generation and analyses of meflA mutant strains. Panel A. Mutant AX2 and AX4 strains were generated by partially deleting the 
mef2A gene through homologous recombination. DNA was isolated from several clones and analyzed by PCR using oligonucleotides specific for 
the mef2A~ deleted region {mefIA) or for ribosomal DNA (rDNA) used as the internal control. The results obtained from non-mutated (AX4, AX2) 
and mutated (clones 37, 2 and 3) samples are shown. Panel B. Wild-type cells (AX4, AX2) and mef2A~ mutant cells {mef2A~), derived from AX4 or 
AX2 cells were clonally grown on K. oerogenes for 4 days. Biological replicates were performed on several different plates. Pictures of the colonies 
were taken, and their size determined in three independent experiments. The average area of the colonies and the standard deviations are 
indicated under each picture. Scale bar: 5 mm. Panel C. Wild-type (AX4, AX2) and mutant cells {AX4/mef2A~, AX2/mef2A~) were collected and 
cultured under starvation conditions to study the multicellular development. The initial steps of aggregation, streaming and mound formation 
were assayed under submerged conditions (Mound column). For later stages of slug (Slug) and fruiting body (Culminant) formation, the cells 
were placed on nitrocellulose filters. Pictures were taken using a Leica stereomicroscope. Scale bar: 0.3 mm. 
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evidence of any transcription of the remaining coding 
DNA of this gene (Gene Expression profile of mef2A 
mutant cells Subsection). The mefIA deleted strains grew 
more slowly than the wild type strains when feeding on 
bacteria (Figure 3B), although no difference in growth was 
observed in the axenic culture (data not shown). 

The mutant strains completed development under 
starvation conditions at the same time as the wild-type 
strains, but several differences were observed during the 
process, as shown in Figure 3C. Both the mutant and 
wild-type strains formed streams during aggregation, but 
the streams of the mutant strain appeared more 
fragmented than those of the wild-type strains, forming 
more heterogeneous and smaller mounds. Subsequently, 
the mutant strains formed fewer slugs, which were 
smaller and migrated shorter distances than those 
formed by the wild-type strains. Finally, the mutant 
strains formed more culminant structures that were 
more heterogeneous in size than those of the wild-type 
strains. The number of spores formed by each strain was 
quantified, and the mutant strains produced about half 
the number of spores formed by the wild-type strains 
(Table 1). However, the viability of the mature spores 
was similar for the mutant and wild-type strains. 

The developmental phenotype of the mutant was 
further characterized by studying the expression of cell- 
type specific marker genes. AX4 and mefIA mutant 
cells were transfected with reporter vectors that drive 
lacZ expression under the control of the ecmA and ecmB 
prestalk gene promoters or the pspA prespore gene pro- 
moter. The lacZ expression was analyzed at the finger, 
slug and mid-culminant stages of development (Figure 4). 
The ecmA gene promoter is active in the anterior, 
prestalk region of finger and slug structures (PstA 
region) and in the stalk, upper and lower cups and basal 
disk of wild-type structures. The activity was detected in 
the same regions in mefZA' mutants, but the anterior 
prestalk region was larger in the fingers. A quantification 
of the relative size of the PstA prestalk region indicated 
that it represented 23.45% (standard deviation: 5.27%) of 
the finger length in AX4 structures and 31.16% (SD 
6.62%) in the mefIA' mutants, difference that is statisti- 
cally significant (p < 0,001 according to the students 
test). In addition, ecmA promoter activity was more 
extended and diffuse in the mutant slugs. ecmB promoter 



Table 1 Production 


of spores by wild type and mef2A' 


mutant strains 








Strain 


Wild 


type 


mef2A' 


AX4 


100± 


13.5 


61.2 ±12.5 


AX2 


100± 


13.9 


41.8 ±15.8 



The value 100 indicates the average number of spores produced in AX2 or 
AX4 cells. The average and standard deviations of three biological replicates, 
made in triplicate, are shown. 



activity in the anterior prestalk region of mound and slug 
structures is more restricted than that of the ecmA pro- 
moter and is mainly located in the tip organizer region, as 
well as in a number of cells scattered around the posterior 
regions, i.e., the anterior-like cells [31]. In culminant struc- 
tures, the ecmB promoter is active in the stalk. The mutant 
strains showed an extended region of ecmB promoter 
activity in the finger structures. In the slug structures, 
ecm^-expressing cells were more disperse in the mutant 
strains, and no well-defined tip organizer region was 
observed. The pspA prespore promoter presented a pattern 
of activity opposite to that of the ecmA promoter. In the 
fingers, the pspA promoter was active in the posterior, 
prespore region. The quantification of the pspA prespore 
region indicated that it is shorter in mef2A' mutants 
(61.32%, SD 9.2%) than in AX4 structures (70.92%, SD 
7.52%) (p < 0,001). In the culminant structures, pspA was 
active in the sorus. The mef2A mutants showed reduced 
pspA promoter activity, especially in the slugs where no 
well-defined prespore region was observed (Figure 4). In 
the culminant structures, the mutant sori were thinner 
than the wild-type sori. 

The differences observed could be due to the participa- 
tion of Mef2A in the process of prespore differentiation 
but could also be due to defective inter-cellular signaling 
in the mutant structures. We designed a developmental 
analysis of mixtures of the mutant and wild-type cells to 
discriminate between these two possibilities. In these 
experiments, we transfected wild-type and mutant cells 
with a reporter vector that would express lacZ upon differ- 
entiation of the cells to prestalk {ecmBdacZ} or prespore 
(pspAdacZ) cells. The transfected cells were mixed with 
non-transfected cells in a 1:4 proportion and allowed to 
develop. If the mefIA' mutant cells were defective in gener- 
ating the intercellular signals required for cell differenti- 
ation, their mixture with wild-type cells would provide the 
defective signal and induce correct differentiation of the 
mutant cells. Alternatively, if the mefIA' mutant cells were 
defective in the process of cell differentiation, the presence 
of wild-type cells would not compensate for their differenti- 
ation defect. Therefore, wild-type and mutant cells express- 
ing lacZ from the prestalk-specific ecmB promoter or the 
prespore-specific pspA promoter were mixed with 
unlabeled cells and allowed to develop. The slugs and early 
culminant structures were analyzed for lacZ expression, the 
results of which are shown in Figure 5. The expected 
distribution of ecm^-expressing cells can be observed in 
the AX^I KXA'lacZ samples. The ecm^-expressing cells 
are located in the tip-organizer region of the slugs, as 
well as scattered in their posterior region. In culminant 
structures, ecmB is expressed in the stalk (including the 
tip) and in the upper and lower- cup regions. The 
mixture of mef2A' mutant labeled and unlabeled cells 
showed a staining pattern similar to that shown in 
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Figure 4 Distribution of prestall< and prespore cells in mef2A mutant developmental structures. Wild-type (AX4) and mef2A mutant cells 
{mef2A~) were transfected with reporter vectors where locZ was expressed under the control of the promoter region from the prestalk-specific 
genes ecmA {ecmAdocZ} and ecmB {ecmBdocZ} or the prespore-specific gene pspA {pspAdocZ}. Pools of transfected cells were placed on 
nitrocellulose filters to induce multicellular development. Structures were collected at the finger, slug and culminant stages of development, and 
/acZ expression was determined by histochemical Xgal staining. Pictures were taken using a Leica stereomicroscope, after counter-staining with 
eosin. Scale bar: 0.2 mm. 



Figure 4 for ecmBdacZ-expressing mutant cells. When 
the ecmB-expressing mef2A cells were mixed with the 
wild-type cells, no lacZ expression was detected in the tip 
organizer region in the slugs or in the tip of the culminant 
structures, indicating that the mefZA cells were excluded 
from these regions. As expected, AX^-ecmBdacZ cells 
were found in these regions when mixed with mef2A' 
cells, indicating that mefZA might be required for prestalk 
cell differentiation at the tip-organizer region. 

The results obtained using pspAdacZ as a cell marker 
are shown in the lower panel of Figure 5. Homogeneous 
mixtures of AX4/AX4 and mef2A' mutant/me/2A' mutant 
cells presented the same pattern of staining shown in 
Figure 4. In the case of the mutant cells, a reduced popu- 
lation of /<3^cZ-expressing cells was also observed in the 
slugs and, to a lesser extent, the culminant structures. 
The mixture of pspA-lsbeled mef2A' mutant cells with 



wild-type unlabeled cells showed that very few mutant 
cells differentiated as prespore cells, and the cells that 
expressed pspA were found at the rear region of the 
slugs and the lower part of the sorus in the culminating 
structures. In contrast, intense lacZ staining was observed 
when pspA-expressing AX4 cells were mixed with mef2A' 
mutant cells. 

Mixing experiments were also employed to study spore 
formation in chimeras. In this case, the cells were labeled 
(Cell-tracker, see MM) and mixed in a 1:1 proportion with 
unlabeled cells. Cells mixtures were set to develop and 
allowed to differentiate for 24 hours. Spores were 
collected, and the percentage of fluorescent spores was 
determined. Table 2 shows that homogeneous mixtures of 
wild-type (AX4 or AX2) or mefZA' mutant cells yielded 
the expected proportion of approximately 50% fluorescent 
spores. However, the mixture of labeled wild-type cells 
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Figure 5 Distribution of mef2A mutant cells, expressing prestalk or prespore markers, in developmental structures formed in 
combination with wild-type cells. Wild-type (AX4) or mefIA mutant {mef2A) cells were mixed in a 4:1 proportion with pools of cells 
transfected with reporter vectors (AX4::/acZ, mef2AdocZ) expressing locZ under the control of the prestalk-specific ecmB promoter (ecmBdocZ) or 
the prespore-specific pspA promoter {pspAwlocZ}. Cell mixtures were allowed to develop on nitrocellulose filters, slugs or culminant structures 
were collected and lacZ expression was determined by Xgal hydrolysis. The two upper rows of pictures show the distribution of cells expressing 
the prestalk ecmB promoter, and the two lower rows show those expressing the pspA prespore promoter. The structures were stained with eosin 
and observed under a stereomicroscope. Scale bar: 0.2 mm. 



with unlabeled mefZA mutant cells produced more than 
90% fluorescent spores. In perfect agreement with this 
result, when labeled mefIA cells were mixed with 
unlabeled wild-type cells, less than 10% of the spores 
showed fluorescence. 

We also analyzed spore formation by in vitro differenti- 
ation. In these experiments, starved cells were induced to 
differentiate into spores by incubation with 8-Br-cAMP 
[32]. After 30 hours, the presence of differentiated spores 



Table 2 Percentage of fluorescent spores in wild-type/ 
mef2A' mutant chimeric structures 





AX4 


AX2 




42.4 ±6.5 


54.9 ± 6.8 


WT^^/mef2A 


91.5± 11.8 


94.4 ±8.9 


mef2A^^/WT 


7.8 ± 2.3 


5.3 ±2.1 


mef2A^^/mef2A 


42.6 ± 7.8 


41.7 ±5.6 



The fluorescent cell population is indicated In the first column by the 
superscript FL for each 1:1 cell mixture. The average and standard deviations 
of 4 biological replicates are shown. 



was determined by the morphological changes observed 
under the microscope (Figure 6). Incubation of AX4 cells 
induced the differentiation of more that 90% of the cells 
into ellipsoid, highly refringent spores. However, when the 
mef2A' mutant cells were treated, the majority of the cells 
appeared rounded and flattened and very few (less than 
5%) refringent spores were observed (lower panel of 
Figure 6). Similar results were obtained for AX2 wild- type 
and mutant cells (data not shown). 

Gene expression profile of meflA mutant cells 

We analyzed the differences in gene expression between 
AX4 and mef2A' strains after 16 hours of development 
(at the finger stage) when the difference in cell-type marker 
expression is greatest (see Figure 4). Poly(A)'^RNAs 
were isolated, converted to cDNA and sequenced using 
an lUumina massive sequencing machine. A total of 
14,847,781 alignable sequences were obtained for the AX4 
RNA and 11,809,110 for the mef2A' RNA, corresponding 
to approximately 13,500 different genes. The sequences 
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Figure 6 In vitro differentiation of wild-type and mef2A mutant 
cells into spores. Wild-type (AX4) and mef2A mutant imef2A~) cells 
were incubated in the presence of Br-cAMP, a cell-permeable 
derivative of cAMP, to induce spore differentiation. Pictures were 
taken after 30 hours of incubation using a TSlOO Eclipse Nikon 
microscope. Scale bar: 10 pm. 

V J 



obtained were aligned over the AX4 genome sequence, 
and the number of sequences obtained for each gene was 
determined. The genes that presented a significant differ- 
ence in the number of reads between wild-type and mefZA 
mutant strains were finally determined. Given that only 
one RNA sample was analyzed for each strain, stringent fil- 
ters were used to determine the genes that were differen- 
tially expressed: more than 3 times in the number of reads, 
with an adjusted p value of less than 0.01. Seventy-seven 
genes showed significant differences with these criteria. 
Thirty two of these 77 genes showed higher expression in 
the wild-type strain, and 45 showed higher expression in 
the mefIA mutant strain. As a control, the meflA tran- 
script was sequenced 73 times in the AX4 sample and 
none in the mutant sample. Table 3 shows the more sig- 
nificant genes that were found coding for known proteins 
or for proteins with a number of conserved domains. The 
largest group of genes code for small proteins that are 
expressed in prestalk cells and that are generally expressed 



to higher levels in the meflA mutant structures. A number 
of these genes showed similarity to the hssA gene [33]. A 
second group of genes coded for 57-59 amino acid long 
proteins that do not show significant similarity to hssA but 
are also expressed in prestalk cells. Three genes coding 
for small proteins (69-72 amino acid long) that were 
expressed in prespore cells showed lower expression in the 
mefIA mutant structures. Other genes that are important 
for prestalk development and that were expressed at higher 
levels in the mep,A mutant structures include Pks32, 
which codes for a polyketide synthase and could be in- 
volved in the synthesis of prestalk differentiation factors, 
and mybCy which codes for a transcription factor involved 
in the response to prestalk differentiation factors [34]. A 
number of genes coding for proteins possibly involved in 
transcription regulation were identified. Their expression 
appears to be dependent on mefIA because the expression 
is significantly decreased in mutant structures. 

To confirm and extend the data obtained, we performed 
quantitative RT-PCR experiments on new mRNAs isolated 
from structures collected after every two hours of develop- 
ment. The expression levels of 11 of the genes listed in 
Table 3 were determined by quantitative RT-PCR, and the 
results are shown in Figure 7. The genes selected were rep- 
resentative of the main categories identified in Table 3 and 
include two /z55A-related genes {}issA, G0283503) and one 
gene coding for a small protein that was expressed at 
higher levels in the wild-type strain (G0285863). The other 
analyzed genes coded for proteins involved in prestalk dif- 
ferentiation {Pks32, mybQ, for proteins possibly involved in 
transcription regulation (G0290847, G0271438, mybQ, for 
spore-inducing factors (psil) and for membrane proteins 
possibly involved in extracellular signaling or cell adhesion 
{tgrFl, tgrCS, G0285697). The results for these genes are in 
complete agreement with the mRNA sequencing data and 
show that the majority of these genes are similarly 
dependent on mefIA for all developmental times analyzed. 
The only exception was psil, which was expressed at later 
developmental stages in the mefIA mutant structures than 
in the wild-type structures (16 hours vs. 6 hours). The ex- 
pression of several of the genes studied (G0285863, 
G0290847, G0271438, G0285697) was almost completely 
dependent on mefIA for all developmental stages analyzed. 

Discussion 

The study of the biological function of the mefIA gene 
was approached by generating deletion mutants in D, 
discoideum. Several mutants were generated in the AX2 
and AX4 strains, and similar phenotypes were observed 
for all strains, as shown in Figure 3 and Tables 1 and 2. 
Over-expression of mefZA caused developmental defects 
that were morphologically similar to those observed in 
the mutants. Pools of meflA over-expressing cells were 
similar to the mutant strains in the formation of more 
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Table 3 Genes that showed a significant difference in their expression between wild-type and mef2A' mutant 
structures developed for 16 hours, as determined by mRNA sequencing 



hssA-related genes (prestalk-specific) 



Small proteins (57-59 aa) (prestalk specific) 



Small proteins (69-72 aa) (prespore specific) 



Tiger family proteins 
Polyketide synthase family 
Transcription regulation 



Developmental genes 
Signaling proteins 



Increased expression in the mutant 

G0267936 -93 (93-1) 
G0268400 -80 (160-2) 
G0277741 -55.6 (1724-31) 
G0281001 -31.2 (780-25) 
G0281 189 -134.75 (539-4) 
G0281191 0 (79-0) 
G0281195 -35.33 (106-3) 
G0281197 -55.5 (111-2) 
G0282307 -106.67 (320-3) 
G0283713 -72.75 (291-4) 
G0293356 -11.24 (281-25) 
hssA. -10 (1550-155) 
G0283421 -34.33 (515-15) 
G0283465 -118.67 (356-3) 
G0283501 -25.19 (3929-156) 
G0283503 -41.42 (994-24) 
G0283505 -33.95 (1935-57) 
G0283507 -38.25 (1224-32) 
G0283511 -9.55 (2015-211) 
G0283515 -14.28 (2313-162) 
G0283519 -55.88 (950-17) 
G02721 88 -15.28 (2119-139) 
G0269674 -15.57 (794-51) 
G0284283 (111-0) 
G0283395 (97-0) 
G0271888 -40.5 (81-2) 
G0269672 -7.46 (574-77) 



tgrF 1.-39.85 (518-13) 
Pks32 -6.62 (1205-182) 
my bC -10.86 (228-21) 



Stl5 -7.69 (1084-171) 
psiN -41.51 (5645-136) 
Omtl2 -19.79 (277-14) 
arrK -10.65 (213-20) 



Decreased expression in the mutant 

G0293362 10.74 (26-279) 
G0281013 11.95 (19-227) 



G0285863 28.07 (208-7523) 
G0284623*c 35.48 (29-1029) 
G0271110* 13.91 (11-153) 
tgrC5*c 18.9 (10-189) 

comH^ca 11.53 (184-2123) 
G0288967* 17.07 (96-1639) 
G0290847*a 21.74 (80-1739) 
G0290855*c 32.39 (41-1328) 
G0271438* 65.5 (4-262) 
srfC. 0 (0-73) 
psir 26.06 (469-12220) 

hspeca 15.69 (26-408) 
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Table 3 Genes that showed a significant difference in their expression between wild-type and mef2A' mutant 
structures developed for 16 hours, as determined by mRNA sequencing (Continued) 



Metabolism 


osbH -38.33 (115-3) 
G0278647 -23.43 (164-7) 


fhbB 6.54 (525-3433) 


Putative Membrane proteins 


G0275535 -6.03 (3197-530) 


G0285697* 451.2 (5-2256) 




G0289143 -9.79 (1116-114) 


G0267564* 16.32 (41-669) 




G0287195 -10.05 (774-77) 


G0272714* 10.99 (73-802) 




G0284683 -11.87 (273-23) 


G0270342* 10.81 (1679-18158) 
G0272042 8.53 (459-3915) 


Translation regulation 




Rpl32* 79.86 (7-559) 
R52*a 10.44(195-2035) 


Other functions 


G0277795 -23.55 (259.11) 


G0284969* 12.87 (93-1197) 
cog2* 12.41 (539-6687) 
G0276325* 32.62 (8-261) 
G0290965*c 12.31 (13-160) 


The following information is indicated in each line of columns 2 and 3: the reference number of the gene, the ratio of the number of reads in wild-type and 
mutant samples (negative when a larger number of reads was obtained for the mutant strain) and the actual number of reads for each sample (mutant-wild type). 
Asterisks on the third column indicate genes that present a sequence related to IVlef2 binding sites in their putative promoter regions. The letter "c" indicates that 
this sequence correspond to the consensus Mef2 binding site. The letter "a" indicates the presence of consensus binding sites for the S. cerevisiae ArgSO 



MADS-box transcription factor. 



Structures, more heterogeneous in size, and fewer slugs 
with reduced migration (data not shown), which prevented 
mutant complementation studies. Similar results were 
obtained when mef2A was expressed from integrative 
vectors using the constitutive ActlS promoter or the 
endogenous, prespore-specific mefZA promoter. A possible 
explanation for these results might be that over-expressed 
MefA proteins bind to co-factors or activating molecules 
outside the chromatin environment, thus impairing their 
regulatory function on DNA-bound Mef2A molecules. 

The results indicate that mefZA is involved in the deter- 
mination or differentiation of prespore cells and of a group 
of prestalk cells in D, discoideum. Mutant cells do not 
differentiate to spores in vitro and, in vivo, produce 
approximately half the number of spores than wild-type 
cells produce. These defects are cell-autonomous because 
the presence of wild-type cells is not able to induce differ- 
entiation of the mef2A mutant cells. The proposed role 
would represent a conserved function for Mef2 proteins 
during evolution, given that plant and animal homologous 
proteins also play important roles in cell differentiation, as 
mentioned in the Introduction. Mef2A mutants also show 
impaired growth when feeding on bacteria. This de- 
fect does not appear to be due to a reduced phago- 
cytic capacity, as determined by incubation with 
fluorescent microspheres (Fluoresbrite, PolySciences, 
data not shown). Differences in cell motility could 
also explain the smaller size of the colonies, but these 
possibilities have not been further studied. However, 
mef2A (srfC) has been previously identified as one of 
the genes whose expression is regulated depending on 



the growth substrate of the D, discoideum cells, bacteria 
or axenic media [35]. 

As mentioned above, mef2A appears to be involved in 
the determination or differentiation of a population of 
prestalk cells located at the tip of the culminant structures 
and at the anterior-most region of the slugs. These cells 
express the ecmB gene and include the tip-organizer cells 
that regulate the culmination of the structures [31]. The 
evidence for this function is that mef2A' mutant cells 
expressing ecmBvlacZ do not participate in the formation 
of the tip in mixed developmental processes (Figure 5). In 
addition, mef2A mutant cells differentiate poorly to ecmB- 
expressing prestalk cells in vitro in the presence of DIF 
and cAMP [36] (data not shown). However, we have not 
detected mef2A expression in these cells. The cell- 
autonomous function of meflA during differentiation may 
be due to a cell-type determination process taking place at 
the previous mound stage of development where mef2A is 
expressed (Figure 2). Alternatively, the expression of 
mef2A in ecm^-expressing prestalk cells might be too low 
to be detected by gene reporter expression analyses. 

Other populations of prestalk cells, however, are enlarged 
in mutant structures, as shown in Figure 4. Differences in 
cell type specification could explain the phenotypes 
observed. For example, the existence of a larger number of 
prestalk cells that adhere more strongly to each other might 
contribute to the breakage of the streams. In addition, a 
number of the genes that are misregulated in mef2A 
mutants, such as the tgr family of genes {tgrFl, tgrCS), are 
involved in cell adhesion [37]. However, no differences 
in Ca-dependent cell adhesion could be determined 
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(See figure on previous page.) 

Figure 7 Developmental expression pattern of representative genes differentially regulated in wild-type and mef2A mutant structures, 
as determined by mRNA sequencing. RNA was isolated from wild-type (AX4) or mef2A~ mutant {mef2A) cells during growth (time 0) or from 
structures developed for the indicated times on nitrocellulose filters (2-20 hours). RNAs were converted into cDNA, and the expression level of 
each gene was determined by quantitative PGR. A fragment of the large mitochondrial ribosomal RNA was used as an internal control for 
expression quantification. The relative value of 1 was assigned to the expression level of growing wild-type cells for each gene. Open bars 
correspond to the wild-type expression levels and black bars to the mef2A~ mutant expression levels. Panels on the left correspond to genes that 
were expressed at the highest levels in mef2A mutant structures, according to the mRNA sequencing analyses, and the panels on the right 
correspond to genes that were expressed at higher levels in wild-type cells. Two biological replicates and three technical replicates of each 
sample were analyzed. Error bars represent the standard deviation of the data. 



experimentally using the method described by Parkinsons 
et al [38] (data not shown). The slug structures also 
showed a highly altered proportion and distribution of 
prestalk and prespore cells (Figure 4), which could explain 
their smaller size and limited motility. 

Despite the developmental defects discussed above, it 
seems clear that mef2A is not absolutely required for D, 
discoideum development given that a large number of 
fruiting bodies and spores are formed in the mutant 
strains. We would like to suggest that meflA participates 
in a network of transcription factors that regulate cell 
differentiation and that could compensate, at least partially, 
for the absence of mefIA, For example, mutations of 
histone deacetylases [39] and chromatin-binding proteins 
[40] affect prespore differentiation or cell-type patterning. 
In addition, the expression of several genes coding for 
putative transcriptional regulators is regulated by mep,A, 
For example, comH codes for a GATA-binding transcrip- 
tion factor expressed in prespore cells. G0288967 codes for 
a putative p-sandwich domain transcription factor, and 
G0290847, G0290855 and G0271438 code for proteins 
containing domains possibly involved in DNA binding and 
transcription regulation. rblA, a retinoblastoma homolog, 
also controls the preference of the cells for stalk or spore 
differentiation [41]. This gene is expressed in the prespore 
region, and, in chimera with wild-type cells, rblA mutant 
cells show a strong preference for stalk differentiation, as 
was also demonstrated here for mefIA mutants. 

Mef2A might play a cell-autonomous regulatory role in 
cell differentiation in response to extracellular signals analo- 
gous to those described in other biological systems. The 
activity of vertebrate Mef2 transcription factors is tightly 
regulated by extracellular signals. One of the best-known 
regulatory pathways involves the regulated association of 
Mef2 with class II histone deacetylases in a process medi- 
ated by protein phosphorylation [42]. Mef2 can also be 
directly phosphorylated through MAP kinase pathways, 
regulating its transcriptional activity [43]. In this respect, 
the D, discoideum MAP kinase ErkB is required for spore 
differentiation [44,45]. It would be of interest to determine 
whether mefIA participates in this ErkB-mediated spore 
differentiation pathway, especially given that Mef2A 
presents three consensus ERK phosphorylation sites. 



The main regulatory pathway described as inducing 
prespore cell differentiation is initiated by extracellular 
cAMP and requires protein kinase A (PKA) activation 
[46]. The in vitro spore differentiation study shown in 
Figure 6 was induced by Br-cAMP treatment, resulting in 
direct PKA activation. The mefZA' mutant cells were 
unable to differentiate under these conditions. This result 
indicates that Mef2A regulation takes place downstream of 
PKA. In fact, a PKA consensus phosphorylation site is 
present close to the C-terminal end of Mef-2, indicating the 
possibility that this protein might be a substrate for PKA. 

Massive-sequencing analyses of gene expression indicate 
that Mef2A can play a role in cell differentiation through 
the regulation of gene expression. There are 32 genes 
whose expression decreases in mefIA mutants, and many 
of these genes are specifically expressed in prespore cells, as 
determined by the mRNA expression analysis available at 
Dictybase (DictyExpress [47]) and by the in situ hy- 
bridization analysis [48]. A number of prespore specific 
proteins, often used as prespore markers, also showed dif- 
ferences in expression between AX4 and the meflA mutant 
strain but did not reach the filter requirements set up in 
the analysis of the sequencing data (more than 3 times the 
difference in the expression level and a p-value smaller than 
0.01). For example, cotA was expressed 2.12 more times in 
AX4 than in the mutant, cotC was expressed 2.02 times, 
cotD 1.48 times, pspD 1.95 times, pspB 252 times and pspA 
1.39 times. In contrast, many of the genes that are over- 
expressed in the meflA mutant have been identified as 
prestalk-specific, as shown in Table 3. The prestalk-specific 
gene ecmA was also expressed at higher levels in the 
mutant (1.51 times). However, ecmB was expressed 1.53 
times more in AX4, in agreement with the reduced and 
more disorganized distribution of ecm^-expressing cells in 
the mutant at 16 hours of development (Figure 4). 

Further studies are required to determine the mechanism 
involved in the transcriptional regulation of the genes 
whose expression is altered in mefIA mutants. Several of 
the genes that are under-expressed in mefIA mutants 
appear to be almost completely dependent on this tran- 
scription factor for their expression. The regulation of the 
expression of these genes could be mediated by the direct 
binding of Mef2A to their regulatory regions. The DNA 
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binding site of Mef2-related factors has been conserved 
through evolution and corresponds to the consensus 
sequence CTA(A/T)4ATG. We looked for the presence of 
this sequence on the 1000 nucleotide-long fragments 
located upstream of the 21 genes down-regulated in the 
mutant using the oPOSSUM program [49]. As shown in 
Table 3, 6 of these genes contained this consensus 
sequence in the region analyzed. Thirteen additional genes 
contained related sequences that differed in the C or the 
G nucleotides. In addition, 4 genes presented the consen- 
sus binding site of a related MADS -box transcription 
factor (ARG80 from S, cerevisiae) [50]. These data indicate 
that a number of these genes could be direct regulatory 
targets of Mef2A. Alternatively, mef2A could regulate the 
expression of other transcription factors controlling the 
expression of these genes. 

The initial analysis of the structure of the mefZA gene 
detected the existence of alternative promoters that drive 
the expression of the gene at different times of develop- 
ment and in distinct structures. It is remarkable that the 
related srfA and srfB genes are also transcribed from 
alternative promoters, specific for different cell types and 
developmental stages [30,51]. Other developmental regula- 
tory genes are also transcribed from alternative promoters 
in D, discoideum, such as pdsA (extracellular phospho- 
diesterase) [52] carA (cAMP receptor) [53] and acaA 
(adenylyl cyclase A) [54]. The existence of alternative 
promoters might have been an evolutionary adaptation that 
regulates the expression of a gene under different condi- 
tions of growth and/or different developmental processes. 

Conclusions 

mefIA, which codes for a protein homologous to myocyte 
enhancer factor 2 transcription factors, is required for 
several of the steps of the D, discoideum biological cycle, 
including growth on bacteria and multicellular develop- 
ment. In particular, mefIA is involved in the regulation of 
the determination or differentiation of prespore cells and a 
group of prestalk cells during the developmental process of 
fruiting body formation. 

Methods 

Cell culture, transformation and development 

D, discoideum cells were cultured axenically in HL5. 
Transformation by electroporation was performed as de- 
scribed previously [55]. Transformed cells were selected by 
treatment with blasticidin [56] or neomycin (G418). Filter 
development was induced by spreading 1-2 x 10^ cells 
(0.6-1.2 X 10^ cells/cm^) on nitrocellulose filters (Millipore 
Co., Bradford, MA, USA) [57]. The phosphate-based PDF 
buffer was used to obtain finger, Mexican-hat and culmin- 
ant structures. Slug structures were obtained by conducting 
the development in water-soaked filters under directional 
light. Submerged development was induced by incubation 



of the cells in 2 ml of PDF phosphate-based buffer 
on 37 -mm diameter cell-culture dishes at 5 x 10^ cells/ml. 

Phylogenetic studies 

The amino acid sequences of MADS -box regions from 
various organisms were obtained from public databases 
and compared to those of the four D, discoideum pro- 
teins containing MADS -box-related sequences, obtained 
from DictyBase (http://www.dictybase.org). Amino acid 
sequences were compared using the ClustalW program 
at the online Biology Workbench facilities from the San 
Diego Supercomputer Center (http://workbench.sdcs.edu) 
and the ClustalX program [58]. Phylogenetic trees were de- 
termined using the neighbor- joining method [59]. A ran- 
dom generator seed of 111 and 1000 bootstrap trials were 
calculated. Trees were drawn using the NJplot program. 

Rapid amplification of cDNA ends 

RNA was isolated from AX4 cells at 8 hours of multi- 
cellular development on nitrocellulose filters. The SMART 
cDNA amplification kit from Clontech (Clontech Laborator- 
ies, Inc, Mountain View, CA, USA) was used for the ampli- 
fication of the 5' untranslated region of mefIA mRNA 
according to the manufacturers instructions. The oligo- 
nucleotide TGTTGCCTGTCTATTTCTTTCATTAG, com- 
plementary to nucleotides 145 to 164 of the gene, was used 
as the primer. Amplification products were cloned in the 
pGEM®-T Easy Vector System (Promega Co, Madison, WI, 
USA), and the insert of at least 10 different colonies of each 
product were sequenced. 

Determination of mef2A expression by RT-PCR 

RNA was isolated from 2 x 10^ cells, either during growth 
or after development on nitrocellulose filters for 2 to 24 
hours, using the TRI reagent (Sigma- Aldrich Inc., St. Louis, 
MO, USA) according to the manufacturers instructions. 
RNA was further purified using the RNeasy Mini Kit 
(Quiagen). cDNAs were generated from 2 [ig of total RNA 
using gene-specific oligonucleotides as primer. cDNAs were 
used as substrates for PGR reactions using as primers the 
oligonucleotides used for cDNA synthesis and upstream 
oligonucleotides designed from the coding region of the 
transcripts. The oligonucleotides GGACTAGTTTCCATT 
GAACCAATTGGGTGAGCG and CTGATAATACAGAT 
AATACTCGC were used for mef2A cDNA synthesis. The 
large mitochondrial rRNA was amplified as a control, using 
oligonucleotides GGGTAGTTTGACTGGGGCGG and 
CACTTTAATGGGTGAACACC. 

Vectors for the generation of knockout strains 

Flanking regions of the mef2A gene, including nucleo- 
tides -4069 to -2692 and 596 to 1298, in relation to the 
A of the translation initiation codon were generated by 
PGR and cloned on both sides of the blasticidin 
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resistance gene in the pLPBLP plasmid vector [60]. The 
ligonucleotides GGCCGCGGCCATTCCCAGCAACGCT 
GGTAATC, GGTCTAGACCTGGAAAACTGGAAAACC 
AATTG, GGATCGATCCACCCACACTAACACACACC 
and GGGTCGACGGTGGTGGTGATTGGTGCTG were 
used for these amplifications. The oligonucleotides TGG 
GAAGGAATAAAATTACAATTGAAAAG and GCGAG 
TATTATCTGTATTATCAG were used to test for the dele- 
tion in AX4 strains, and GTTGCCTGTCTATTTCTTTC 
and CACTCACTTACATATCACACACC were used in 
AX2 strains. 

Construction of the reporter and expression vectors 

The two mef2A promoter regions were amplified by PGR 
from D. discoideum genomic DNA and cloned in the 
reporter vector PsA-ialphaGal [61] in substitution of the 
Xbal/Bglll pspA promoter fragment. The oligonucleotides 
GGTCTAGAGCACAAGATTATACTTGCCA and GGA 
GATCTCATGGTGTGTGATATGTAAGTGAGTG were 
used to amplify the -2201 to -489 region of the gene, 
corresponding to Promoter 1. The oligonucleotides 
GGTCTAGACACTCACTTACATATCACACACC and 
GGAGATCTTGTTGCCTGTCTATTTCTTTCATTAG 
were used to amplify the -511 to 164 region, correspond- 
ing to Promoter 2. The complete promoter region (-2201 
to 164) was amplified using the first and last oligonucleo- 
tides described above. Previously described lacZ reporter 
vectors were used to determine the expression of the 
developmental markers pspA [61], ecmA and ecmB [62]. 

The mef2A gene was expressed using the pDV- 
CGFP-CTAP vector [63], under control of the ActinlS 
promoter. The region coding amino acids 4 to 1046 
of the protein, including the third intronic region of 
the gene, was amplified from genomic DNA using the 
Pfx DNA Polymerase (Invitrogene ). Two overlapping 
fragments were obtained using the oligonucleotide 
pairs GGATCCAGGAATAAAATTACAATTGAAAAG/ 
GGAGATTGATGCTGTGGTTG and CAACAACAAAG 
CGCCAATCC/ACTAGTAGGTTCCATTGATTTTCTTT 
TTCGG. The two fragments were joined together 
using the overlapping Haell restriction site and the 
resultant fragment cloned between the BamHI and 
Spel restriction sites of the vector. A second expres- 
sion vector where mef2A was expressed under control 
of his own promoter was constructed by substituting 
the ActinlS promoter of the pDV-CGFP-CTAP/me/2A 
vector by the mef2A promoter. The ActinlS promoter 
was excised by Sail and BamHI digestion and 
replaced by the complete mef2A promoter previously 
cloned in the PsA-ialphaGal vector, as described 
above. The mefZA promoter was isolated by Xbal/ 
Bglll digestion. The Sail end of the vector and the 
Xbal end of the promoter were converted to blunt 
ends before ligation. 



Histochemistry and determination of p-galactosidase 
activity in developmental structures 

Cells transformed with the different reporter vectors were 
allowed to develop on nitrocellulose filters for the time 
periods indicated in each experiment. The structures 
were fixed and permeabilized, and |3-galactosidase ac- 
tivity was detected by hydrolysis of X-gal (S-Bromo-4- 
chloro-3-indolyl |3-D-galactopyranoside), as previously 
described [64]. 

Cell tracking experiments 

Growing cells were collected by centrifugation and 
resuspended in phosphate-based PDF buffer containing 
S (il/ml of 100 mM CellTracker™ Blue CMHC (4- 
chloromethyl-7 -hydroxy coumarin) (Invitrogen, Eugene, 
Oregon, USA) or a vehicle (DMSO) and incubated for 1 
hour in shaking cultures. Cells were then washed, 
resuspended in free PDF buffer and mixed in a 1:1 propor- 
tion. A total of 6.6 X 10^ cells from each mixture were 
spread out on nitrocellulose filters for 24-36 hours, after 
which several sori from each mixture were harvested and 
dissociated in water. Spores were visualized in a Zeiss 
Axiophot fluorescence microscope and counted. 

Determination of mRNA levels by quantitative RT-PCR 

RNA was isolated from 2 x 10^ cells, either during growth 
or after development on nitrocellulose filters for the times 
indicated in each experiment, using the TRI reagent 
(Sigma- Aldrich Inc., St. Louis, MO, USA) according to the 
manufacturers instructions. RNA was further purified 
using the RNeasy Mini kit (Quiagen). cDNAs were gener- 
ated from 2 (ig of total purified RNA using random primers 
(Promega Co., Madison, WI, USA). cDNAs were used as 
substrate for quantitative real-time PGR reactions using 
the following gene-specific oligonucleotides: hssA gene 
(DDB_G0280999) GTGCTATTACCTCAATTTCAAG and 
GGCAACCACATGAACCACTTG; DDB_G0283S03 gene 
CAAATCATTACAATCAATCACAAGTG and GGGCTA 
CAGCAGCAACTG; prSl gene (DDB_G028S863) CCAA 
TAATTCTTTGAAGGCCC and CAATAGCTTGGCCC 
ATAGTAGC; tgrFl gene (DDB_G0292732) CCCACCA 
TTTACTCCAATACTC and GTAGAGATGGTGTTGAT 
GGAG; tgrCS gene (DDB_G0281407) GCTGGCTTAGC 
ACTTTCATCAG and GAGACCAACGGCAGCGACAC; 
pks32 gene (DDB_G0292732) CAACTCCAGTCACAA 
CTATAGC and GATTATCATGAATGTGGAATGCTG; 
mybC gene (DDB_G0281S63) GGTGGAGGTAAAACTG 
GTGC and CATCCATCCAACTAATATCACG; DDB_G 
0290847 gene CAGTACTGAACAAGCATTATCAAG and 
GTTAACATAACCTTGTTGAGAATC; DDB_G0271438 
gene GTCATGAAATTGGAGATCGAAG and CATGAGA 
TGATGTTGATTTGG; psil gene (DDB_G0288919) GGT 
TGTACACTTGTACCACG and GAGGTGCTTCAAAG 
AGAGC; DDB_G028S697 gene GGTAAGGCAGTTGTC 
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N: Polyasparagine tract; Q: Polyglutamine tract. Parcel B. The amino acid 
sequer^ces of the MADS-box domair^s (amino acids 1 to 60 of the SrfC 
sequence) and the contiguous SRF- or Mef2-specific domains (amino acids 
61 to 86) of the species indicated on the left are aligned in relation to the D. 
discoideum SrfC (l\/lef2A) sequence using the Clustal W program. Darker 
boxes indicate the presence of 8 or more identical amino acids. Lighter 
boxes indicate the existence of 8 or more amino acids with similar chemo- 
physical characteristics. Figure S2. Analysis of the structure of the 5' region 
of the mef2A gene and determination of the transcription start sites. The 
sequence of the 5' region of the mef2A mRNA was determined by primer 
extension using the rapid amplification of the cDNA ends (RACE) technique. 
The sequences obtained were aligned with the genomic DNA sequence to 
determine the intron/exon structure of this region of the gene and the 
transcription initiation sites [GenBank:KC852901]. Exon regions are indicated 
in black capital letters while intron sequences are indicated in blue small 
letters. Consensus splicing sites are underlined. The sequence is numbered 
from the Adenine of the translation initiation codon, shown in bold 
characters. Transcription initiation sites, as determined by RACE, are 
indicated with red asterisks over the nucleotide sequence. The complete 
sequence of the third intron of the gene and the sequence of the large 
fourth exon, coding for the C-terminal region of the protein, are not shown 
but are schematically indicated on Figure 2. 



AATGC and GCCTACCAGCTGAGACTTCAGC. A re- 
gion of the large mitochondrial ribosomal RNA was ampli- 
fied as a loading control using the oligonucleotides 
CACTTTAATGGGTGAACACC (used as a reverse oligo- 
nucleotide) and GGGTAGTTTGACTGGGGCGG (used as 
a forward oligonucleotide). The StepOnePlus Real-Time 
PGR System (Life Technologies Co., Applied Biosystems, 
Carlsbad, CA, USA) was used in these experiments. PGR 
products were labeled with SYBR Green using the Power 
SYBR® Green PGR Master Mix reaction mix (Applied 
Biosystems) following the manufacturers instructions. The 
final volume of the reaction was 20 (il, using a 0.2 [iM con- 
centration of each primer. PGR conditions were as follows: 
95°G, 10 m; (95°C, 15 s; 45°G, 30 s; 62° 1 m) x 30-40 cycles. 

In vitro spore differentiation 

Exponentially growing cells were washed in KK2 buffer 
(16.5 mM KH2PO4, 3.9 mM K2HPO4, 2 mM MgS04, 
pH 6.2), plated on culture dishes at a concentration 
of 10^ cells/ml in spore buffer (10 mM MOPS, 20 mM 
KGl, 20 mM NaCl, 1 mM CaG^, 1 mM MgCl2, pH 6.2) 
and supplemented with 12.5 mM 8-Br-cyclic-AMP and 
20 \iM GdCl2 [32,65]. Gells were incubated in the dark 
for 30 hours and then observed under a TSlOO Eclipse 
Nikon microscope (Nikon, Tokyo, Japan). Pictures were 
taken with a Leica DFG420 camera (Leica Microsystems, 
Wetzlar, Germany). 

mRNA sequencing 

RNA was isolated from structures developed on nitrocellu- 
lose filters for 16 hours using the TriReagent and purified 
with an RNeasy Mini kit, as described previously. Poly 
(A) -containing RNA was isolated and converted to 
cDNA. The cDNA was fragmented, amplified by PGR 
and the nucleotide sequences determined using an 
lUumina Genome Analyzer IIx massive sequencer at 
the Parque Gientifico de Madrid. Sequencing data 
were analyzed at Sistemas Genomicos, S.L. (Valencia, 
Spain). The generated sequences were mapped to the 
D, discoideum genome using the TopHat vl.1.3 software 
[66]. Transcripts were identified and quantified using 
the Gufflinks vl.0.3 program. The total number of 
reads per gene was determined using the HTSeq 
package (http://www-huber.embl.de). Statistical analyses 
of the results was performed using the DESeq package 
[67], using an FDR of 0.01. A minimal difference of three 
times in expression levels was considered. 

Additional file 



Additional file 1: Figure SI. Comparative study of the amino acid 
sequer^ces of the D. discoideum MADS-boxcontair^ing proteins. Panel A. The 
functional and structural domains present in the D. discoideum SrfC (Mef2A) 
protein are schematically shown. M; MADS-box; m: Mef2-conserved domain; 
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